A Fitter Code for Deep Virtual Compton Scattering 
and Generalized Parton Distributions. 
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We have developed a fitting code based on tlie leading-twist handbag Deep Virtual Compton Scat- 
tering (DVCS) amplitude in order to extract Generalized Parton Distribution (GPD) information 
from DVCS observables in the valence region. In a first stage, with simulations and pseudo-data, we 
show that the full GPD information can be recovered from experimental data if enough observables 
are measured. If only some of these observables are measured, valuable information can still be 
extracted, with certain observables being particularly sensitive to certain GPDs. In a second stage, 
we make a practical application of this code to the recent DVCS Jefferson Lab Hall A data from 
which we can extract numerical constraints for the two H GPD Compton form factors. 

PACS numbers: 



I. MOTIVATION 

Generalized Parton Distributions (GPDs) have 
emerged during the past decade as a powerful concept 
and tool to study nucleon structure. They describe, 
among many other aspects, the (correlated) spatial and 
momentum distributions of the quarks in the nucleon 
(including polarisation information), its quark- antiquark 
content, a way to access the orbital momentum of the 
quarks, etc. 

Formally, in short, the GPDs are Fourier transforms of 
matrix elements in Quantum Chromo-Dynamics (QCD) 
for light-cone bilocal operators between nucleon states of 
different momenta. For helicity conserving quantities in 
the quark sector, there are four GPDs, H, H, E, E which 
depend, in leading order and leading twist QCD, upon 
three variables: x, ^ and t. Both x and ^ express the lon- 
gitudinal momentum fractions of the two quarks of the 
bilocal operator, while t is the squared four-momentum 
transfer between the final and initial nucleon. Experi- 
mentally, GPDs are most simply accessed through the 
measurement of exclusive leptoproduction of a photon 
(Deep Virtual Compton Scattering -DVCS-) and, possi- 
bly, of a meson. GPDs arc then accessed through the 
factorization with an elementary perturbative process, 
which leads to the so-called handbag diagram, predicted 
to be dominant at small t and large Q^, where Q^ is the 
virtuality of the initial photon. Fig. [1] shows schemati- 
cally these notions for the DVCS process on the proton 
on which, as a first step, we will concentrate our work 
and discussion in this article. Indeed the leading twist 
contribution is expected to be the most directly accessi- 
ble for DVCS and proton targets currently provide (and 
will also provide in the near future) the richest data set. 
We refer the reader to refs. [i|,[3,[3, 4, 5, 6, 7] for the origi- 
nal theoretical articles and recent comprehensive reviews 
on GPDs for more details on the theoretical formalism of 
GPDs and DVCS. 

Concerning actual data, the field is now rapidly grow- 
ing. The first data related to DVCS on the proton were 
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FIG. 1: The handbag diagram for the DVCS process on the 
proton ep —> e'p'y' . Here x + ^ and x — ^ are the longitudinal 
momentum fractions of the initial and final quark, respec- 
tively, and t = {p ~ p')^ is the squared momentum transfer 
between the initial and final protons (or equivalently between 
the two photons). There is also a crossed diagram which is 
not shown here. 



the beam spin asymmetries (BSA) which have been pub- 
lished by the HERMES 8] and CLAS ^ collaborations 
in 2001. Although the kinematic range over which the 
data were integrated was large and, to some extent, the 
actual exclusivity of the reaction could be disputed, these 
were the first very encouraging signals that the handbag 
diagram could be accessed, since several theoretical pre- 
dictions were in relatively good agreement with the data. 
Then, longitudinally [l3| and transversely ll] polarized 
target s pin asymmetries as well as beam charge asym- 
metries |12| were measured which, though suffering from 
the same experimental limitations as the BSAs, also con- 
firmed the general theoretical expectations. 



With these first encouraging signals found in non- 
dedicated experiments, a second generation of DVCS ex- 
periments has been launched in the past few years, with 
fully dedicated detectors, i.e. high resolution electromag- 
netic calorimeters at Jefferson Lab (JLab) to detect the 
final-state photon and a recoil detector at HERMES to 
detect the recoil proton. This has resulted in the pub- 
lication of the first proton DVCS cross sections (beam- 
polarized and unpolarized) in the valence region at a few 
precise kinematical points by the JLab Hall A collabora- 
tion [1J| and several BSAs over a large kinematic range 
by the JLab HaU B collaboration ^. HERMES results 
with the recoil detector are now eagerly awaited. 

The field looks extremely promising with a series of 
novel experiments planned, in particular at JLab, aimed 
at accurately measuring new DVCS observables, includ- 
ing: 

• Longitudinally polarized target spin asymmetries 
and cross sections (along with double polarizations 
observables) with JLab experiment E05114 ^15i |. 

• Transversely polarized target spin asymmetries and 
cross sections (along with double polarisations ob- 
servables) with JLab experiment E08021 [iq . 

• New precise (unpolarized and beam-polarized) 
cross sections at new kinematics with JLab exper- 
iments E06003 [3 and E07007 ^. 

A similar experimental program in the longer term fu- 
ture is planned with the JLab 12 GeV upgrade, which 
will cover a larger phase space (see JLab approved ex- 
periments E1206114 [11] and E1206119 [20|]). 

Anticipating this rich harvest of data, one asks the 
question how to extract the GPD information from them. 
We recall that the DVCS process is accompanied by the 
Bethe-Heitler (BH) process, in which the final state pho- 
ton is radiated by the incoming or scattered electron and 
not by the nucleon itself. The BH process, which is not 
sensitive to GPDs, is experimentally indistinguishable 
from the DVCS and interferes with it, thus complicat- 
ing the matter. Considering that the nucleon form fac- 
tors are well-known at small t, the BH process is however 
precisely calculable theoretically. 

Another issue is that the GPDs, which are a function of 
the three variables x, ^ and t, enter the DVCS amplitude 
within a convolution integral over x. Therefore, only ^ 
and t are accessible experimentally (in the Bjorken limit, 
^ ~ i-a /2 ^"^ which xb is the standard Bjorken vari- 
able). Formally, the DVCS amplitude is proportional to: 
/-I '^^ a;±f=rt'e '^ •'■ (''^licre thc ellipsis stands for similar 
terms for E, H and E). Decomposing this expression into 
a real and an imaginary part, we find that the maximum 
information that can be extracted from the experimental 
data at a given (^,t) point is H{±^,^,t), when measur- 
ing an observable sensitive to the imaginary part of the 
DVCS amplitude, and J_, dx — f l , when measuring 



an observable sensitive to the real part of the DVCS am- 
plitude. 

If ones reduces the range of x from { — 1, 1} to {0, 1} in 
the convolutions, there are in principle eight GPD-rclated 
quantities that can be extracted: 
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We will call these, respectively, in a symbolic notation, 
Re{H), Re{E), Re{H), Re{E), Im{H), Im{E), Im{H) 
and Im{E). These are also often called the Compton 
Form Factors (CFFs). Note here the absence of — tt fac- 
tors in our definition of the Im{) CFFs with respect to 
ref. [21|. In the following, for convenience and in a very 
loose way, we will speak of real part CFFs and imagi- 
nary part CFFs to designate generically the CFFs de- 
fined by eqs. [l]|l]and by eqs. [S][H1 respectively and which 
correspond to the CFFs associated with the real and 
thc imaginary parts of the DVCS amplitudes, respec- 
tively. The CFFs can be decomposed into terms for 
individual quark flavors which, for the proton, yields: 
H{U,t) = |i/"(C,e,0 + ^H^i^^tt) and similarly for 
E, H and E. The quark flavor separation, which wc will 
not tackle in this study, can be carried out by measur- 
ing DVCS on a neutron target, which yields a different 
quark flavor combination. The first measurements of the 
neutron BSA have been recently published |22l | . 

In summary, given the interference of the BH process, 
the deconvolution issue regarding x and the large number 
(8) of independent quantities to be extracted from the 
data, it is clearly a non-trivial task to extract the GPDs 
from the experimental data and, ultimately to map them 
in the three variables x, ^, t. 

The first stage of any general program of measuring 
GPDs should certainly be to extract the eight CFFs from 
the data for a given ^,t point, in a model-independent 
way. This would be only the beginning of the program, 
since the x dependence would still need to be deconvo- 
luted using in principle a model with adjustable parame- 
ters. At this stage, let us mention that there might actu- 
ally be a couple of ways around this issue: firstly, if the 



Double-DVCS process, (i.e. with a virtual photon in the 
final state) can be measured, then varying the virtuality 
of the final state photon provides an extra lever arm and 
this allows us to measure the GPDs at each x,(_,t val- 
ues (thoueh with some limitations if the final photon is 
timelike) [23|, |2j| . Secondly, dispersion relations could in 
principle reduce from eight to five the number of GPD 
quantities to be extracted, by expressing the real part 
CFFs defined by eqs. [l]|4] in terms of integrals over ^ of 
the respective imaginary part CFFs defined by eqs. [Mil 
plus a real subtraction constant (at fixed ^ and t). This 
strategy requires one to measure data over a very wide 
range in ^ (at fixed t) unless one has good reasons to trun- 
cate the inte gral or to extrapolate. We refer the reader to 
refs. [25I . [26l . \2T\ . [28| for discussions on this subject and, 
in particular, to ref. [23| for actual fits to the DVCS data 
in collider kinematics. 

As a first approach, in this article, we focus our study 
on the most general case, in which the eight quantities 
of eq. [T][8] are independent. Our present purpose is to 
understand to what extent, given the leading twist and 
leading order QCD DVCS amplitude, the CFFs can be 
extracted from various observables. We present our work 
in three steps, corresponding to the next three sections. 
In the next section, we present the general framework 
of the fitting code that we have developed. In particu- 
lar, we present the tests that we have carried out in a 
Monte-Carlo approach using simulated data. With ideal 
(i.e. insignificant) statistical error bars, we will learn 
general features such as which observable is sensitive to 
which CFF, what are the highest reconstruction efficien- 
cies that one can expect to achieve in such ideal condi- 
tions, etc. In the following section, we will simulate real 
experimental conditions by introducing a smearing of the 
simulated data together with realistic error bars. We will 
then discuss the resulting reconstruction efficiencies and 
the uncertainties on the reconstructed GPD parameters. 
Finally, in a last stage, in section [TVl wc will apply our 
fitting code to the recent JLab Hall A DVCS data and 
will attempt to extract quantitatively some of the first 
real GPD quantities. 



II. MONTE-CARLO STUDY 

In this section, we present the general features of our 
fitter code and, in particular, the way we have tested it 
and established its reliability and efficiency. Our fitter 
code is simply based on the merging of the well estab- 
lished VGG [23, ISO, [Slj code, which calculates numeri- 
cally, using model CFFs, the leading order and leading 
twist handbag DVCS + BH amplitudes and observables, 
and the well known MINUIT minimization program from 
CERN (33. The VGG code uses its own models for the 
GPDs but in fact any GPD model can be used. For our 
fitting purposes, the idea is to consider the CFFs that 
enter the DVCS amplitude as free parameters to be fit- 
ted and to see how well they can be extracted from the 



DVCS Monte-Carlo (MC) data. 

We have started by testing and studying our code with 
pseudo-data. The procedure we have followed is based on 
a MC approach and can be summarized by the following 
general steps; for a given experimental kinematic point 
uniquely defined hy E^, Q^, ^ and t [E^ being the beam 
energy): 

1. We generate randomly a set of values for the 
quantities Re{H), Re{E), Re{H), Re{E), Im{H), 
Im{E), Im{H) and Im{E), 

2. From these values, we calculate various observables 
which we will detail shortly as a function of 0, the 
azimuthal angle between the leptonic and hadronic 
planes, 

3. We fit with our code these distributions with eight 
"parameters" , which are meant to correspond to 
Re{H), ReiE), Re{H), Re{E), Im{H), Im{E), 
Im{H) and Im{E), 

4. Wc finally compare the fitted values to the gener- 
ated ones and draw our conclusions. 

This procedure is going to be looped over several hun- 
dred times, exploring the whole phase space of values, 
within given limits, that can be taken by Re{H), Re{E), 
Re{H), Re{E), Im{H), Im{E), Im{H) and Im{E). In 
doing so, we do not make any assumption or do not take 
into account any information on the absolute or the rela- 
tive values of the generated GPDs, such as, for example, 
whether Im{H) is dominant over Im{H) and Im{E), 
which could well be the case in reality. In this way, we 
place ourselves, in the framework of this very general 
study, in the "bfindest", the least biased and the most 
conservative of situations. Entering physically motivated 
information or contraints on the CFFs can of course only 
improve the results that we obtain. Before presenting 
the results of our Monte-Carlo study, let us mention a 
few details. 

We will consider nine independent observables which 
can be expected to be measured with relatively good ac- 
curacy in the near future. These are: cr, Acr^O: Acto^, 
Ao-Qy, Actoz, Acr2a:, Aa^y, Aazz and Aac- Here a refers 
to the unpolarized cross section. When there are two 
indices. Act refers to the difference of polarized cross sec- 
tions and the two indices refer respectively to the polar- 
ization of the beam and of the target (i.e. Act^q is the 
difference of cross sections with a longitudinally polarized 
beam and an unpolarized target, as has been measured 
recently by the JLab Hall A collaboration JJJ] , and AcTzz 
corresponds to the difference of cross sections that can be 
obtained with a longitudinally polarized beam and a lon- 
gitudinally polarized target). Finally, Actc refers to the 
difference of unpolarized cross sections between a nega- 
tive and a positive lepton beam, i.e. proportional to the 
beam charge asymmetry measured by HERMES |I3 |. 

Regarding the parameters to be fitted, as a first ap- 
proach, we have taken into account only seven CFFs, 



instead of eight: Re{H), Re{E), Re{H), Re{E), Im{H), 
Im{E), Im{H), i.e. we have set Im{E) to 0. The reason 
is that the GPD E is usually associated with the pion 
pole i-channel exchange which is real. Nothing keeps us 
in principle from considering Im{E) and taking it as an 
extra free parameter in the fit. However, there is also 
clearly no need to complicate the fit procedure if it is 
not justified. Conversely, we have decided to let Re{E) 
be a free parameter, even though it is supposed to re- 
flect the pion pole and is therefore well-constrained. In 
this way, the conjecture that E comes from the pion pole 
can be verified. If this conjecture proves false, the whole 
parametrization of E will clearly have to be revisited in- 
cluding its imaginary contribution. Thus Im{E) = 
is the only model assumption that we will make in this 
study. 

Also, we have decided to take as fit parameters, not 
the CFFs of eq. [l][8] themselves, but the deviations with 
respect to them. This is not an assumption, simply a con- 
vention. In other words, we take some (supposedly) real- 
istic reference values for Re{H), Re{E), Re{H), Re{E), 
Im{H), Im{E) and Im{H) and fit the coefficients that 
multiply these reference values to the MC data. There- 
fore in our MC study, we are going to generate ran- 
domly, in a given range, seven real numbers that mul- 
tiply the seven reference CFFs. From these new (ran- 
dom) CFFs, we generate the DVCS observables. We fit 
them and our aim is to recover the seven initial mul- 
tiplicative coefficients, knowing the numerical values of 
the reference CFFs. We will call these multiplicative 
coefficients, the GPD multipliers and denote them as 
a{Re{H)), a{Re{E)), a{Re{H)), a{Re{E)), a{Im{H)), 
a{Im{E)) and a{Im{H)). Similarly to the CFFs, in the 
following, for convenience, we will speak of real part and 
imaginary part multipliers, although, mathematically, all 
these numbers are obviously real. These can be inter- 
preted as "ratios" of the fitted CFFs to the reference 
CFFs. Clearly, this is only a matter of convention and 
instead of these GPD multipliers, we could certainly have 
generated and fitted Re{H), Re{E), etc. directly. Our 
general motivation is to create the most efficient and ro- 
bust code by starting the minimization procedure as close 
as possible to the true solution. This implies that the ref- 
erence CFFs must be as realistic as possible. For these 
reference CFFs, we have taken, as a first approach, those 
given by the VGG parametrization. It is by no means 
infered that the VGG parametrization is the most real- 
istic one that currently can be found in the literature, 
but this was obviously the most convenient choice in our 
case because they provide a decent description of several 
existing DVCS experimental data @, [E[3- However, it 
is clear that instead of VGG as a reference, any other 
GPD model could be taken. 

In the present study, we have generated seven numbers 
a{Re{H)), a{Re{E)), a{Re{H)), a{Re{E)), a{Im{H)), 
a{Im{E)) and a{Im{H)). The range we have considered 
is {—4,4}. This means that we allow the randomly gen- 



erated CFFs to vary in absolute value up to a factor 4 
from the VGG value. This range {—4,4} is quite arbi- 
trary and docs not have any significant impact on the 
subsequent results of this section. Whether such large 
variations from the VGG CFFs are realistic or not is not 
particularly the object of this section which is meant as 
an exercice, testing the ability to extract GPDs. Let us 
just mention that the GPDs of the VGG code are normal- 
ized and respect the various model-independent relations 
{Hi{x, 0, 0) = q{x), H^ix, 0, 0) = Aq{x), where q{x) and 
Ag(x) are the unpolarized and polarized parton distri- 
butions functions of flavor q, respectively) and sum rules 
{j^ldxH'i{x,^,t) = F^{t), etc.). Although these rela- 
tions are quite constraining and do not permit a fully 
arbitrary normalisation of the GPDs, let us note that 
there are contributions to the GPDs which can escape 
any formal normalization constraints. For instance, the 
so-called D-term [33], which cancels at ^=0 and is odd in 
X so that it does not contribute to the form factor sum 
rule, is such a contribution. Also, GPD parametrizations 
based on other approaches such as in ref. [3J| can find in 
the valence region an H GPD larger by a factor 2 to 3 
compared to VGG. 

Once the seven initial multipliers are generated, the 
next step is to calculate with the VGG code the 
DVCS-I-BH amplitude and the various (unpolarized, 
singly polarized, doubly polarized) resulting cross sec- 
tions. We recall that, in this study, we restrict our- 
self to leading- twist. However, in order to ensure gauge 
invariance of the DVCS amplitude, we mention that 
the VGG amplitude contains a small (on the order 
of 10% depending on the kinematics) twist-3 contribu- 
tion [29.. .30.] . Since Wandura-Wilczek kinematical twist- 
3 [35| and some other contributions to higher twists 
are also available in the VGG code {k± effects for in- 
stance) [23, 130| , let us mention that fits can also be done 
beyond leading twist at the cost of longer computing 
time. 

At given E^, ^, Q^ and t values, we generate the cj) dis- 
tributions of the nine observables previously mentioned. 
Fig. m shows such distributions calculated from the set of 
randomly generated a() multipliers which are displayed 
in table [D The kinematics correspond to that of the 
recent JLab Hall A data [ll|: Ee=5.75 GeV, a;s=0.36 
(corresponding to C= 0.22), Q2=2.3 GeV^ and -t=0.28 
GeV^. Our goal is therefore to recover the a() values 
of table U from the fit to the distributions of fig. [2l 
Once again, we recall that these values are the coeffi- 
cients that muhiply the VGG CFFs, i.e. not the CFFs 
themselves. The corresponding VGG CFFs are also dis- 
played in table HI The VGG code has several options for 
the parametrization of the GPDs. As we said, it does 
not really matter which particular option is used since 
the VGG CFFs are just meant to be used as starting or 
reference values and what is ultimately extracted are the 
deviations from these values. We have used the Regge 
inspired unfactorized ansatz for the i-dependence '311 of 
H and E, a factorized t dependence ansatz for H and the 
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FIG. 2: The nine independent DVCS observables considered 
in this work as a function of (j>, at the kinematics Ee=5.75 
GeV, a;B=0.36, Q^=2.3 GeV^ and -f=0.28 GeV^ The curves 
have been generated with the GPD muhipliers of table HI The 
left column, from top to bottom, corresponds to the unpolar- 
ized cross section with an electron beam (solid curve) and a 
positron beam (dashed curve) and the three target single po- 
larization differences of cross sections: Actqz, Aaox and Aaoy. 
The right column, from top to bottom, corresponds to the 
beam polarized difference of cross section Aa^o and the three 
beam-target double polarization differences of cross sections: 
Aa,,, Aa,x and Aa, 



All observables are four-differential 
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pion pole for the modelling of E. Also, no D-teiuv was 
used and the parameters h^jai = hsea = 1 in the profile 
function of the Double Distributions were selected. 

Obviously, since the GPD multipliers of table |T] have 
been randomly generated for our study, they have no 
particular meaning, and the curves of fig. [5] have no re- 
lation to the actual a and Actzo that the JLab Hall A 
collaboration has measured at similar kinematics. 

In this section, we generate the (j) distributions with 
"ideal" error bars, i.e. 5x10""' in relative value to the 
cross sections. The idea is first to gain confidence in our 
fitting program as well as to understand, under ideal con- 
ditions, some general features of the fit procedure such as 
which GPD is sensitive to which observable, the number 
of observables needed to extract all GPDs, and so forth. 
In the following section, we will discuss real conditions 
by assigning realistic error bars to these distributions. 

In fitting these 4> distributions with the seven free 
parameters a{Re{H)), a{Re{E)), a{Re{H)), a{Re{E)), 



a{Im{H)), a{Im{E)) and a{Im{H)), we stress that we 
ignore any information that we know about the initial 
generated values, the idea being to match as much as 
possible real conditions. In this study, we set the start- 
ing values to 0. For real data, as we said, in order to 
be as close as possible to the true solution, the starting 
values should be set to 1. The range over which the val- 
ues of the seven coefficents are allowed to vary is set to 
{—5, 5}, i.e. somewhat larger than the range of the gen- 
erated values. This point is essentially the only piece of 
information that we take from the initial input. We have 
used a least square method and the MIGRAD minimizer 
of MINUIT to perform the fit. The quantity that we 
minimize is thus: 
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where atteo is the theoretical cross section (or difference 
of cross sections) from the VGG code, aexp is the corre- 
sponding (simulated) experimental value (generated from 
the VGG code as well) and 6aexp is its associated ex- 
perimental error bar. We recall that in this section, we 
consider ideal error bars. The index i runs over the num- 
ber of observables to be fitted, i.e. nine at maximum in 
the present context. For all the results that we present in 
the next two sections of this article, fits have been carried 
out with seven experimental points in </> spread in steps 
of (/)=30° between 0° and 180°. Little improvement is ob- 
served in fitting more </> values however, computing time 
increases. Fitting seven experimental points for each of 
the nine observables with our seven parameter formula- 
tion takes about thrity minutes, running on an average 
computer. Fitting only a and Actzq takes less than ten 
minutes. 

We have generated several hundred events in order to 
have statistically significant results where an event is 
understood to be a set of seven real values a{Re{H)), 
a{Re{E)), a{ReiH)), a{Re{E)), a{Im{H)), a{Im{E)) 
and a{Im,{H)). Let us now see how the fitted values 
of the GPD multipliers compare to the generated ones. 
We present our results with the support of two figures. 
In the following, we call a topology a particular combi- 
nation of some of the nine observables: <t, Atr^O: Actok, 
AtToj/, AtToz, Actzj;, Acr^y, Act^^ and Actc- Fig. [3] com- 
pares the generated and the reconstructed values of our 
seven GPD multipliers under different topologies, using 
one particular event corresponding to the generated GPD 
multipliers of table HI Fig. [5] summarizes our results for 
several hundred events, in the phase space {—4, 4} of the 
a() multipliers, for each of the fitted parameters, under 
the different topologies of fig. [31 The particular event of 
fig. [3] has been chosen for illustrative purposes because it 
allows us to make several (not all, though) of the conclu- 
sions that we can draw with hundreds of events. Let us 
start by discussing this figure. 

In fig. [31 the solid lines show the (randomly) gener- 
ated value of the seven multipliers (see table [l|, while all 





ReiH) 


Re(E) 


Re{H) 


Re{E) 


Im{H) 


Im{E) 


Im{H) 


a() 

VGG value 


0.378 
0.658 


-1.818 
0.355 


3.296 
0.458 


-0.699 
41.705 


-3.732 
1.58 


2.608 
0.48 


-1.285 
0.43 



TABLE I: One set of seven randomly generated GPD multipliers a() which are used in the illustration of our study, together 
with the corresponding VGG values for the seven CFFs at a;s=0.36 (corresponding to (,— 0.22), and — i=0.28 GeV^. 



the points show the fitted (reconstucted) values obtained 
from our fitting program, under different topologies. We 
are going to discuss them one by one, in parallel with 
figure m to understand if the features that we can infer 
from fig. [3] are general or simply an accident. 
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FIG. 3: One particular example of the fit results for different 
topologies. The solid line indicate the generated values of the 
a() GPD muhipliers of table U 

In fig. [21 the empty circles show the reconstructed val- 
ues obtained if only one observable, the DVCS BSA -° , 
(which is so far the most widely measured observable) is 
fitted. We see in this figure that only a{Im{H)) is well 
reconstructed, as well as a{Im{H)), to a somewhat lesser 
extent. On very general grounds, it is not at all surpris- 
ing that fitting only one (/> distribution, i.e. the BSA 
here, with seven parameters, we are not able to recover 
all of them. From fig.lH we see that the good reconstruc- 
tion of a{Ira{H)) is a relatively general feature while the 
reconstruction of a{Im{H)) is more of a statistical acci- 
dent. Fig. m shows, for several hundred randomly gener- 
ated events, the ratio of the number of well- reconstructed 







FIG. 4: Reconstruction efBciencies for all the topologies of 
the seven GPD multipliers and of all of seven in concert. The 
symbols coresponding to the topologies have been defined in 
fig. |3l Right bottom plot: reconstruction of all GPD multi- 
pliers simultaneously. The reconstruction efficiency is defined 
as the ratio of the number of well reconstructed events to the 
number of generated ones. An event is "well reconstructed" 
when the reconstructed numerical value matches the gener- 
ated one to within 5% (dotted curve), 15% (dashed curve) 
or 25% (solid curve). The dash-dotted curve for the Im{H) 
BSA corresponds to the reconstruction efficiency (at the 15% 
level) if the starting value for Im{H) in MINUIT is set to the 
generated one ±20%. 



events to the number of generated ones, which we call a 
reconstruction efficiency, for each GPD multiplier. "Well 
reconstructed" means that the reconstructed numerical 
value of a multiplier falls within 5% (dotted curve), 15% 
(dashed curve) or 25% (solid curve) of the generated 
value. One learns from this figure that by fitting only 
the BSA one recovers, at the 25% level, about 55% of 
the generated values of a{Im{H)) while only « 30% and 
less than 20% if the reconstructed values are within 15% 



and 5% of the generated ones, respectively (first column 
of the upper right plot of fig. [4|) . Once again, the efficien- 
cies for recovering a{Im{H)) are calculated with all a()'s 
in the range {—4,4}. Therefore, ii a{Im{H)) happens to 
be close to zero while the other a()'s not, it shouldn't be 
surprising that a{Im{H)) is not well-recovered. If in real- 
ity Im{H) dominates over the other GPDs then these ef- 
ficiencies can only increase. Therefore, the numbers that 
we quote for these efficiencies and which are displayed 
in fig. m should be taken with great care because they 
really depend on the relative value of each CFF. In addi- 
tion, the kinematics can also change the weight of each 
CFF. Thus these efficiencies reflect more a "sensitivity" 
of a given CFF to some observable{s) rather than an ab- 
solute efficiency and they have only a general statistical 
meaning in the present context. 

Another way to improve the efficiency of the fit is to 
set the starting values close to the true solution. It will 
then be less likely that MINUIT falls into some local 
minimum. For instance, the dash-dotted curve in the first 
column (i.e the BSA) of the upper right plot (i.e Im{H)) 
of fig. [4] shows the reconstruction efficiency (at the 15% 
level) if the starting value in MINUIT for Im{H) is set to 
the generated Im{H) (smeared at the 20% level) instead 
of at as we have done so far. Biasing and guiding the 
fit in such a way more than doubles the probability of 
recovering the original value. 

Despite all these caveats, we can still learn a lot from 
fig.m We can note that a{Im{H)) is not recovered more 
than 15% of the time, which means that the good re- 
construction noted in fig. [3] is not generally the case and 
that the reconstruction of the Im{H) CFF in this topol- 
ogy is thus very dependent on its particular value and 
weight with respect to the other CFFs. On the con- 
trary, fig. m shows a feature that fig. [3] might not have 
hinted to, i.e. that the BSA has some non-negligible 
sensitivity to Re{H) since in « 35% of the cases, it is 
well-reconstructed at the 25% level. 

The sensitivity of the BSA to Im{H) is a well-known 
feature since it can be shown analytically (see ref. [2l| for 
instance) that Im{H) is a dominant contribution to the 
difference of beam polarized cross sections Actzo which 
is the numerator of the BSA. The (limited) sensitivity of 
the BSA to Re{H) can be understood through its con- 
tribution to the denominator of the BSA which is the 
unpolarized cross section a where the real part CFFs en- 
ter. 

In fig. [21 the solid circles show the reconstructed values 
which are obtained for our seven parameters, if now the 
two observables just discussed, the unpolarized cross sec- 
tion a and the difference of beam polarized cross section 
AcTzOi are both used in the fit rather than just their ra- 
tio. In fig. [21 we see similar features to the previous case 
where only the BSA was fitted, i.e. that only a{Im{H)) 
is fairly well recovered: -3.354 reconstructed against - 
3.732 generated (table |T|. It also confirms that the good 
reconstruction of a{Im{H)) in the case of the BSA was 
more an accident than a solid feature. However, the ad- 



vantage of fitting simultaneously these two observables 
instead of simply their ratio is shown well in fig. [H where 
we see that the reconstruction efficiency oi a{Im{H)) now 
reaches more than 85% (at the 15% precision level) while 
it was less than 35% in the BSA topology. This illus- 
trates the strong correlation between Acr^o and Im(H) 
and that Im{H) can be recovered, through this topol- 
ogy essentially independently of the values of the other 
GPDs. This conclusion that Acr^o is mostly sensitive to 
Im{H) and barely to the other GPDs is not new and was 
already pointed out for instance in ref. 2l| analytically. 
We confirm this numerically. 

The topology marked with an asterisk in fig.[31is the re- 
sult of the fit when adding one more observable. In addi- 
tion to the unpolarized cross section a and the difference 
of beam polarized cross section Act^qj we now also fit the 
difference of longitudinally polarized target cross sections 
Aaoz- One sees in fig. [31 that the four real part CFFs are 
as poorly reconstructed as in the previous cases. How- 
ever, now there is a second imaginary part CFF, besides 
a{Im{H)), that is well-recovered: a{Im(H))=-l.Q09 in 
the fit compared to the generated value -1.285 (table [l|. 
We see from fig. [H that this conclusion, i.e. the strong 
sensitivity of Actoz to Im{H), which was also already 
pointed out in ref. [2l|, is very general. In more than 
80% of the cases, a{Im{H)) is well reconstructed in this 
topology. We also note that the precision on a{Im{H)) 
is significantly improved. The reconstruction efficiencies 
corresponding to the 5% (dotted curve) and 15% (dashed 
curve) criteria have increased (up to « 80% compared to 
w 50% in the previous topology for the dashed curve). 

We now bring two more observables into the fit, i.e. 
the difference of cross sections obtained with the trans- 
verse X and y transverse polarizations of the target. The 
corresponding results are shown by the open diamonds 
in fig. [31 We observe that the three imaginary part CFFs 
are now well-reconstructed and that this conclusion re- 
mains true generally. Fig. [51 shows that in more than 
80% of the cases, a{Im{H)), a{Im{H)) and a{Im{E)) 
are well-reconstructed in this topology. We can thus in- 
fer that the target transverse (single) polarisation observ- 
ables are very sensitive to Im{E) and that measuring the 
five observables a, Aazo, Auoa:, Acroy, and Actoz allows us 
to extract the three imaginary part CFFs reliably. Fur- 
thermore, fig. m shows that, although it is not illustrated 
particularly well by fig. [31 the combination of these five 
observables has a significant sensitivity to the four real 
part multipliers. About 40% of them are reconstructed 
well. 

Next, the upright solid triangles show the results of fit- 
ting only the beam charge difference of cross sections (i.e. 
no polarization observables). We see now that the three 
real part multipliers a{Re{H)), a{Re{E)), a{Re{H)) are 
reconstructed well while no imaginary part multiplier is 
recovered at all. Fig. [4l confirms that this is the case 
at basically the 100% level (the few well-reconstructed 
imaginary part multipliers in this topology in fig. [31 are 
clearly accidental). Because, on general grounds, Actc is 



expected to be sensitive only to the real part of the DVCS 
amplitude, this result comes as no surprise. We also note 
in fig.[3]that the real part multiplier a{Re{E)) is not re- 
constructed in this topology. This is a general feature 
since, as seen from fig. HI a{Re{E)) is not reconstructed 
in any of the cases with our high statistics event sample. 
Again, the couple of percent of events well-reconstructed 
clearly result by coincidence. This feature can be demon- 
strated analytically by the same kind of argument used to 
show that the difference of beam polarized cross sections 
Acr^o is insensitive to Im(E). 

In fig. [21 the upside-down solid triangles show the ef- 
fect of fitting Acr^o and Actoz, in addition to Actc- Now, 
in addition to the three real part multipliers a{Re{H)), 
a{Re{E)), a{Re{H)), the three imaginary part multipli- 
ers are well reconstructed following our previous discus- 
sions. We note that there is a slight decrease of the re- 
construction efficiency of the three real part multipliers 
when fitting simultaneously Auc, Act^q and Actqz com- 
pared to fitting only Actc. This is due to the fact that, 
here, we have minimized the sum of the three x^'s cor- 
responding to the minimization of Actc, Act^q and Actqz, 
and it can be shown that the derivatives with respect 
to the imaginary part CFFs are the highest compared 
to those with respect to the real part CFFs. In other 
words, in doing so, one loses some sensitivity to the real 
part CFFs. Therefore, in such a situation, where a set 
of observables is sensitive only to some other indepen- 
dent CFFs, like Actc to the real part CFFs, and another 
set of observables only to some other independent CFFs, 
like AcTzo and Acoz to the imaginary part CFFs, instead 
of fitting simultaneously both sets, one should obviously 
adopt a two-step procedure and fit independently each 
set of observables. To illustrate this trivial effect, we 
leave in fig. 3] the results of our simultaneous fit, with 
the understanding that in this topology, if we had fitted 
separately Actc, we would have recovered the same ef- 
ficiencies of the upright solid triangles topology for the 
real part CFFs. This being said, we then see that six of 
the seven GPD multipliers can be reconstructed in a very 
economical way by measuring only three observables, and 
in particular, that the gold-plated way to access the real 
part CFFs is with the beam charge difference of cross 
sections. 

However, the next symbol (solid square) in fig.[3]shows 
that there is another way of accessing the real part mul- 
tipliers with high efficiency. The solid-square topology 
contains the double-polarization observable Acr^x in ad- 
dition to the open-diamond topology, which was made up 
of the unpolarized cross section and the four single beam 
or target polarization observables. We then see that all 
real part multipliers (in particular including a{Re{E)) , to 
which Auc was insensitive) and imaginary part multipli- 
ers are well recovered. Fig. |3] confirms that the particular 
event of fig. [3] is not an accident and that these conclu- 
sions can be generalized. In about 80% of the cases, the 
imaginary part multipliers are well-reconstructed at the 
25% level, similarly for higher precisions, and in more 



than 65% of the cases for the real part multipliers. 

The solid stars show the results obtained when fit- 
ting all polarization observables simultaneously, i.e. now 
adding the two remaining double polarization observables 
A(Tzx and Act^j, to the previous topology. We see from 
fig.[3]that all seven parameters are well-reconstructed, as 
in the previous topology, but fig.[4]shows that this is now 
at a higher efhciency which reaches « 90% for all seven 
multipliers. Fitting eight observables with seven parame- 
ters overconstraints the fit, with redundant information, 
and therefore significantly improves the efhcicncy. 

The remaining symbols in fig. [3] will be discussed in 
the next section. So far, we have described the gen- 
eral features of our fitting program in the ideal situa- 
tion, i.e. with perfectly precise data, and have demon- 
strated its general reliability and power to recover all 
of our seven arbitrary randomly generated parameters 
(CFFs) from fits to the (p distributions, given enough ob- 
servables. We learned which observables are sensitive to 
particular CFFs and that, if only some of these observ- 
ables are available, we could still extract some specific 
CFF. In particular, one can access in principle the real 
part CFFs with relatively high efficiency using double 
polarisation observables in addition to beam charge dif- 
ferences. 

We now turn to realistic conditions. Ultimately, this 
program has to be used on real data with finite preci- 
sion. In the next section, we therefore simulate real ex- 
perimental conditions and find how robust and reliable 
the program remains. 



III. A MORE REALISTIC STUDY 

We are now going to simulate realistic pseudo-data by 
assigning experimental error bars to the data to be fit. 
The procedure consists not only in assigning an error bar 
to the simulated (f) data points but also in smearing the 
central value according to a gaussian probability distri- 
bution whose standard deviation is equal to the error 
bar. 

Inspired from refs. [IJ, [lM|, we will assign, at all (/)'s, 
a 3.5% error bar to the unpolarized cross section. For 
all the differences of cross sections, we will consider two 
cases: 10% or 5%, the former being approximatively the 
error bar on Acr^o which is quoted in refs. [IJ, [la] and 
the latter being a slightly more challenging experimental 
goal. 

We return to fig. [3l The open square and open star 
symbols show the result of the two last topologies dis- 
cussed in section [m (i.e. the solid square and solid stars) 
when realistic error bars are applied to the pseudo-data 
in the distributions. In each of the seven plots, there 
are two open squares and two open stars. The leftmost, 
for both symbols, is the result of the fits when the differ- 
ences of cross sections are smeared by 10% and the right- 
most when they are smeared by 5%. This can clearly be 
inferred from the size of the error bar of the fitted pa- 



rameters which is systematically smaller when the data 
are smeared by 5% compared to 10%. The uncertain- 
ties on the fitted parameters that are presented here are 
the quadratic errors from MINUIT, which is sufficient in 
the present context where we are looking at statistical 
effects. However, we will see in the next section when 
fitting the JLab Hall A a and Acr^o real data that the 
error determination requires a more dedicated study, es- 
pecially when dealing with a underdetermined problem 
with more parameters to fit than available observables. 

We see in fig. El for this particular event, that the gen- 
eral conclusions that we reached in the previous section 
remain, i.e. that essentially all seven GPD multipliers 
are relatively well-recovered when fitting the six observ- 
ables (T, Acr^o, AcToa:, Actoi/, Actoz and Actzz (open square 
topology) and eight observables with IS.Uzx and Act^j, in 
addition (open star topology). However, while in the 
previous section, the original central values were recov- 
ered almost perfectly in most of the cases, there is now a 
clear dispersion of the reconstructed central values with 
respect to the generated central ones. We note that this 
difference is always within two standard deviations and 
therefore that the errors estimation provided by MINUIT 
seems very reasonable. 

Fig. [S] confirms this conclusion with a larger sample 
(though limited in the display in order not to overcrowd 
the figure) of events where we observe that for basically 
all open stars events (with 10% error bar for the pseudo- 
data) the difference between the generated and recon- 
structed values Aa() is within three ct^jq's, with a^i\ be- 
ing the MINUIT uncertainty previously discussed. The 
few cases where it is not the case can easily be identified 
after inspection; either the x^ of the corresponding fit is 
very bad or the generated value is close to 0, meaning 
that the observables had very little sensitivity to it. 

Fig. [5] also shows the solid circles resulting from fits of 
only a and Actzq- We note that these symbols are visi- 
ble essentially only in the Re{H) and Im{H) plots. This 
confirms the observations that we made in the previous 
section that this topology is sensitive to only these two 
GPD quantities (predominantly Im{H)). Indeed, the ab- 
sence of solid circles on the other plots means that their 
uncertainty (TaQ is off-scale, indicating the total lack of 
sensitivity of this topology to these GPD quantities. For 
this topology, we also observe the larger UaO values in av- 
erage for Re{H) than for Im{H) giving evidence for the 
stronger sensitivity of the solid circle topology to Im{H). 

For the open star topology, the very small <T(^a) val- 
ues that one can observe for the a{Im{H)), a{Im{H)), 
a{Re{E)) multipliers, compared to the a{Re{H)), 
a{Re{E)), a{Re{H)) and a{Im{E)) multipliers reflect 
their very high reconstruction efficiency and reliability. 
This is confirmed by fig. [5] which shows for the open- 
star topology, the reconstruction efficiencies of each GPD 
multiplier. The three GPD multipliers just outlined have 
clearly the highest reconstruction efficiencies. Compar- 
ing the open stars of fig. [5] to the solid stars of fig. IH 
we note a loss of reconstruction efficiency when realistic 



simulated data are fitted, which is particularly important 
for Im{E) and Re{E). For these two GPD multipliers, 
the reconstruction efficiency is smaller by a factor of 2 or 
more. 
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FIG. 5: Uncertainty on the fitted parameters cr^o, as given 
by MINUIT, as a function of the difference between the fitted 
and generated values, for the open-star (with 10% error bar 
for the pseudo-data), solid-circle and open-circle topologies. 
The dashed line indicates the three o'(a) limits. 



In Fig. [5l the open circle symbols (i.e. the fit to only 
—2^) show up essentially only in the Re{H) and Im{H) 
plots. However, the difference with the solid-circle topol- 
ogy (i.e. the fit to a and Atr^o separately) is that most 
of the points are now NOT within the three cjaQ bands, 
in spite of the relatively small CToQ values. This can lead 
to dangerous interpretations where we might believe that 
the fitted value of, say Im{H), are extracted with high 
precision because of the small aaQ while the value is not 
necessarily reliable. Firstly, this suggests that a detailed 
and dedicated error analysis must be done when dealing 
with underdetermined systems such as this one where one 
fits only one observable by seven parameters. We will de- 
tail this in the next section. Secondly, let us emphasize 
that this in no way means that BSAs are dangerous or 
not useful, but that either they should be complemented 
by some other observables that provide additional con- 
straints or that some physics input (such as the domi- 
nance of H over the other GPDs) must guide the fitting 
procedure. We recall that we have blinded ourselves, in 
this study, by not telling the fit which values it should 
start with. 



Let us also note in fig.[3]that the error bars of the open 
star topology are always smaller than the ones from the 
open squares, which is consistent with the fact that fit- 
ting eight observables rather than six improves the qual- 
ity of the fit. Making a blunt generalization from this 
particular event one can also conclude that it is almost 
equivalent to fit eight observables with a precision of 10% 
(left, open stars) as to fit six observables with a precision 
of 5% (right, open squares). 
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FIG. 6: Reconstruction efficiencies of three "experimental" 
topologies for the seven GPD multipliers. Left: efficiency of 
the fit of a and Aa^o with the experimental cp distributions 
being smeared by 3.5% for a and 10% for Aa^o- Center: 
efficiency of the fit of a, Aa^o, Aaox, Aaoy, Aaoz, Aa^x, 
Affzy and Aa^z with the experimental <f) distributions be- 
ing smeared by 3.5% for a and 10% for all the differences 
of cross sections. Right: same as center but with the dif- 
ferences of cross sections being smeared by 5%. The dotted 
histograms correspond to a reconstruction efficiency within 
5% , the dashed ones within 15% and the solid ones within 
25%. 

In this section, we have simulated real experimental 
conditions by (Gaussian-)smearing, with a standard de- 
viation equal to realistic experimental error bars, the ob- 
servables that we fit with our code. We conclude that 
the fitting procedure is, though somewhat less efhcicnt, 
still very reliable (the E CFFs suffering the most impor- 
tant loss of reconstruction efficiency) and that the general 
conclusions that we drew in section |TT] remain valid. Hav- 
ing gained confidence in our code and fitting procedure 
through these simulations with pseudo-data, we now pro- 
ceed in the next section by applying the code to some real 



data. 



IV. APPLICATION 

The JLab Hall A collaboration has recently measured 
a and Acr^o UM at a;B=0.36 and Q^=2.3 GeV^ for four 
different t values. From the previous sections, we have 
learned that the simultaneous fit of these two observables 
allows us to access Im{H) at a relatively high efficiency 
level and to a somewhat lesser extent a few other CFFs 
such as Re{H). 

We have therefore run our fitting code using the JLab 
Hall A data. Fig. [7] shows our results. The numerical 
results of the fit for each GPD multiplier are displayed 
in tabic |TT1 

Several issues are to be discussed concerning this table. 
Let us start with the uncertainties on the fitted param- 
eters. In the previous section, the uncertaintites quoted 
were the standard quadratic errors from MINUIT, i.e. es- 
timated from the second derivative of the x^ function to 
be minimized with respect to the parameter under con- 
sideration, based on the assumption that the x^ function 
is parabolic near its minimum. This was sufficient for 
our general purpose of statistically estimating global re- 
construction efficiencies and general features of our code. 
Now, in the present case with real data, the determi- 
nation of the uncertainty on the fitted parameters has 
to be refined. Indeed, we are facing a multidimensional 
problem (seven parameters to fit) which 1) is non-linear, 
2) has potentially strong correlations between the pa- 
rameters and 3) is severely underconstrained with only 
two observables to fit seven CFFs. Therefore, a care- 
ful and detailed error analysis was carried out with the 
MINOS subpackage of MINUIT which allows us to ex- 
plore in a gradual and automated way the x^ landscape 
around the minimum and define one-standard deviation 
uncertainties for each parameter when it reaches x^ + 1- 
This method, which is essential for non-linear problems, 
yields asymmetric error bars. These are presented with 
the symbols a~,. and a\ (the negative and positive er- 
rors, respectively) in table HIl The price one pays is that 
MINOS is very time-consuming and that the error bar de- 
termination of our seven parameters took several hours 
of computing time. 

In table |TT1 the symbol cx3 means that the uncertainty 
coming out of MINOS is very large, i.e. exceeding the 
limits given in MINOS. To be precise, the results in ta- 
ble [ni have been obtained by setting, like in the previ- 
ous section, the lower and upper limits of the domain of 
variation allowed for the seven parameters to -5 and 5 
respectively. The oo symbol means that x^ -I- 1 was never 
reached in the interval considered. It can be interpreted 
as a whole range of values for the associated parameter 
can accomodate the fit with relatively equally good x^- 
In other words, the corresponding GPD multiplier (or 
equivalently CFF) is essentially unconstrained and there- 
fore no particular confidence and meaning can be given 



E =5.75 GeV, x„=0.36, Q^=2.3 GeV^ 




(^ (deg.) 



<ti (dcg.) 



<ti (dcg-) 



<tl (dcg.) 



FIG. 7; Top row: DVCS a, and bottom row: DVCS Actzo as a function of <^, at 2;_b=0.36 and Q^—2.'i GeV^ for four different t 
values. The data are from the JLab Hall A collaboration 13|. The curves are the results of the fit from our code to the a and 
AiTzo observables. The fitted GPD multipliers which produce these curves are displayed in table HTl 
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TABLE II: Fitted GPD multipliers a{) and their negative (c^o) and positive (o"^(0 uncertainties resulting from the fit of the 

JLab Hall A a and Actzo observables. The "x^ /Ndof value for these fits is presented in the rightmost column. As explained in 
the text, only central values with finite negative and positive uncertainties have a reliable meaning. 



to the quoted numerical value. This does not mean that 
the values associated with an cx) uncertainty in table |TT] 
can take any value, irrespective of the other multipliers. 
Many of the seven parameters are highly correlated with 
other parameters and changing the value of one parame- 
ter will change the values of others. This means that for 



those parameters with an oo uncertainty, the solution is 
not unique and that other sets of values are also possible. 
For instance, fixing some GPD multipliers to particular 
values 0, or 1 (the associated CFF takes its VGG value 
in this latter case) and letting all other parameters re- 
main free can result in fits with an almost equivalent x^ 



compared to letting all seven parameters free. 

However, what is remarkable is that even though many 
CFFs have an cx) uncertainty in table HH a few of them, 
depending on their t values, come out with finite error 
bars. For instance, a central value with finite negative 
and positive uncertainties can be extracted for a{Im{H)) 
at all t values and for a{Re{H)) at the largest three t val- 
ues. For these multipliers, the uncertainties, while finite, 
can be very large, the negative errors ranging from more 
than 100% at the smallest t values to « 20% at the largest 
t values. This should not come as a surprise, given the 
complexity and the underdetermination of our problem. 
Therefore, the large uncertainties that we obtain don't 
reflect a lack of quality or precision of the experimental 
data but rather a lack of constraints and of sensitivity 
of the observables to certain parameters. In these con- 
ditions, it is already a significant success to be able to 
pull out some reliable and stable information from such 
partial inputs. We restate that in this study we have left 
all seven GPD multipliers (or equivalently CFFs) free. If 
some of them can be fixed, neglected or constrained to 
a domain smaller than {—5,5}, it is clear that these un- 
certainties can drastically be reduced. This will be the 
subject of another forthcoming article. We recall that 
when enough experimental observables are available to 
fit, as we have shown in the previous sections with sim- 
ulations, such assumptions won't be necessary and all 
parameters should be uniquely determined. 

We have carried out many checks to get confidence 
in the numbers (central values and error bars) associated 
with finite uncertainties in table|TTl For instance, to start 
the minimization, MINUIT requires us to set some start- 
ing values for the parameters to be fitted. As we men- 
tioned in the previous sections, our strategy is in general 
to set them to the VGG values in order to be, hope- 
fully, close to the true solution. However, we checked 
that by starting with randomly generated starting val- 
ues, MINUIT -l-MINOS were always converging towards 
the same central values given in table [ll] and producing 
the same corresponding uncertainties at the few percent 
level. 

Another check was to vary the limits of the domain 
of variation of the parameters (initially set to {—5,5}). 
We made the fits with the domains {—3, 3} and {—7, 7}. 
The outcome was that, although all central values asso- 
ciated with oo uncertainties in table HIl could be different, 
the few values associated with finite positive and nega- 
tive uncertainties (i.e. a{Im{H)) and a{Re{H))) were 
found remarkably stable at the few percent level. It was 
interesting to note that for the two intermediate t values 
a{Re{E)) and a{Re{H)) would systematically take the 
lower limit of the domain (i.e. -3 or -7, similarly to the 
-5 in tablcini though this had no effect on the central val- 
ues oi a{Im{H)) and a{Re{H)). However, a (moderate) 
effect was found for the uncertainties of a{Im{H)) and 
a{Re{H)). It was observed that their uncertainties, cal- 
culated by MINOS, were slowly increasing as the limits 
of the domain of variation increased. This can certainly 



be attributed to the fact that the correlations between 
the "stable" a{Im{H)) and a{Re{H)) multipliers on the 
one hand and the "unstable" a{Re{E)) and a{Re{H)) are 
not completely absent, and that larger (absolute) values 
of the latter multipliers (since they actually reach the 
limits of the domain of variation as mentionned earlier) 
can certainly tend to increase the error bars of the former 
ones. A few paragraphs below, we make a quantitative 
comparison of this (moderate) effect. 

Although the uncertainties on the a{Im{H)) and 
a{Re{H)) are rather large, we can notice a few trends. 
a{Im{H)) is always consistent with 1 within error bars 
and it tends to increase with t. Although the purpose of 
this study is not at all to prove (or disprove) the VGG 
model, it shows that the VGG model seems to provide a 
reasonable estimate of Im(H), possibly underestimated 
at the « 20% level. The parameter a{Re{H)) has a 
strong tendancy to rise with t and although the uncer- 
tainties are large, it seems to be significantly different 
from the VGG prediction and it can differ by up to a 
factor of 4. 
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FIG. 8: Upper plots: the a{Im{H)) (left) and 
a(i?e(7y)) (right) GPD multipliers as a function off extracted 
from the fits of the JLab Hall A data (fig.[7| at ^= 0.22. Lower 
plots: the resulting Im{H)) (left) and i?e(ii") (right) OFF val- 
ues. For a{Im{H)) and Im{H), three results for each t are 
presented corresponding to fits with the domains of variations 
{—3,3}, {—5,5} and {—7,7} from left to right respectively. 
The dotted boxes around the largest t results indicate that the 
associated fit has a x^ about a factor 2 worse (see table |TI]) 
than the fits at the other t values. 
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TABLE III: The reference VGG CFFs at ^= 0.22 for the four JLab Hall A t values. The muhiplication of these reference VGG 
values by the fitted GPD multipliers of table |ll] (and their associated error) yield the measured CFFs. 



We display in table lllll the VGG reference values which, 
by multiplying the a() GPD multipliers and their associ- 
ated errors of table [TTl allow us to extract the measured 
values of the CFFs themselves. We plot in fig.[8]the mea- 
surements of the two multipliers a{Im{H)) and a(Re{H)) 
and the associated CFFs for the t values for which they 
have finite uncertainties. In this figure, we show three 
results for each t for a{Im{H)) and Im{H). These cor- 
respond to the fits with the domains of variations {—3, 3}, 
{—5,5} and {—7,7} from left to right, respectively. The 
stability of the central value as well as the moderate in- 
crease of the error bar with the increase of the limit of 
the domain previously mentioned can be observed. As 
deviations of CFF values by a factor 7 (or even 5) from 
VGG predictions are not easy to conceive, the uncertain- 
titles quoted in table |IT] and displayed in fig. [5] are very 
conservative. They could be reduced if, based on some 
models or educated guesses, smaller limits could be set 
on the other CFFs. 

One should also note from table [TTl that all fits have 
a x^ close to 1 except for the bin at t = —0.33 GeV^. 
We remark that a problem with this same bin was also 
observed in ref. [36| where, in an analysis mainly de- 
voted to the JLab Hall B BSAs, the JLab Hall A BSAs, 
i.e. "''" , were fitted by the phenomcnological function 
1 f^'" .X where the B coefficient was found to take a pe- 
culiar value for this particular t bin, in clear distinction 
from the surrounding t bins. Also, in the present anal- 
ysis, although the results for a{Re{E)) and a{Re{H)) 
are not very meaningful because of their oo uncertain- 
ties, one can nevertheless clearly observe a discontinuity 
in their fitted values as they jump rather abruptly from 
a largely negative values at the three low t values to a 
largely positive value for the last t. Therefore, there def- 
initely seems to be a change of behavior in the data at 
this large t bin. A bad x^ may be indicating that our 
theoretical modeling is not appropriate. We recall that 
our analysis has been carried out at the twist-2 handbag 
level and that this intriguing change might signal the rise 
of higher twist corrections. 

Let us also mention that still a few more pieces of in- 
formation can be extracted from the values of table [TTl 
For several t values, a{Re{E)) has finite uncertainties. 
Since these error bars are systematically over 100%, it 
makes the information difhcult to use, but we can sim- 
ply note that the value of 1, which corresponds to the 
VGG pion pole, is never excluded and that our original 



assumption to neglect Im{E) on this basis is therefore 
not ruled out. There are a few cases in table [TTl where 
only one uncertainty {(t~,-. in general) is finite, which still 

yields some information. In particular, Im{H) for all t's 
has a reliable (lower) bound. 

We have also tried to fit only the Acr^p observables (i.e. 
without the simultaneous fit of a) which would be sensi- 
tive only to Im{H), Im{E) and Im{H). In general, the 
fitted central values for these CFFs were unstable, unless 
we guided and constrained the fit by setting limits on the 
domain of variation of the parameters. In general, we 
obtained values of Im{H) about 30% higher than those 
quoted in table [Hi but they corresponded to much larger 
values of Im{H), compared to those of table [TTl which 
with its opposite sign would provide some compensation. 
This kind of unstability is not unexpected because fit- 
ting only one observable with three parameters is even 
less constraining than fitting two observables with seven 
parameters. 

Finally, we tried to fit the numerous JLab Hall B 
BSAs [3| • We found that without any guidance of the fit 
or any physics input, no fit stability and robustness could 
be reached. This could actually be anticipated from our 
studies in the previous section. However, as we saw with 
fig. [21 guiding the fit by setting the initial values in MI- 
NUIT to values close to the expected ones, which can be 
relatively safely done for Im{H) following our results on 
the Hall A data that we just presented, could tremen- 
dously enhance the fit efficiency and stability. We delay 
to a forthcoming publication the results of such guided 
or educated fits. 



V. CONCLUSION 

In sumary, we have developed a fitter program to ex- 
tract GPD information from various DVCS observables 
at QCD leading twist and leading order. We have gained 
confidence in the program by checking its reliability and 
efficiency on simulated, ideal or realistic, pseudo-data. 
We came to the conclusion that by fitting enough un- 
polarized, singly and doubly polarized observables, with 
realistic experimental error bars, the code was able to 
extract in most cases the seven CFFs considered in this 
work, even in a very blind framework where all param- 
eters are left free and essentially unconstrained. Apply- 



ing realistic and educated constrains, such as dispersion 
relations or model motivated ansatzes, can only reduce 
the number of independent parameters or limit the range 
of variation of certain parameters and therefore improve 
the efficiency and reliability of this fitter program. We 
recall that most of the aforementionned observables will 
be available in the near future from various experiments 
at JLab in particular. If only some of these nine observ- 
ables are available, we have demonstrated that valuable 
partial GPD information could still be extracted: in par- 
ticular, with (T, Aa^o and Actq^, which are planned to 
be measured in the very near future, the Im{H), Im{E) 
and Im{H) CFFs could be very reliably be known. For 
the real part CFFs, one can either measure the beam 
charge difference of cross sections or double polarization 
observables. 

Finally, we have extracted from the JLab Hall A 
data first numerical values with associated uncertain- 
ties for the H GPD combinations H{£_, S,, t) - H{-£_, S,, t) 

and P ^l dx[H{x,S,,t) ~ H{-x,^,t)]C+{x,C)- The cor- 
responding error bars are rather large because we are 



fitting more parameters than observables and no a priori 
knowledge of any CFFs has been assumed in this study. 
Ultimately, when enough observables are available, this 
should indeed not be needed, and we can consider that 
we have placed ourselves in this long-term perspective. 
In a forthcoming article, we will discuss strategies to re- 
duce these uncertainties on the present data if educated 
assumptions are made about the CFFs. This code is 
only the first step towards a general fitting procedure 
of DVCS data (and potentially, if data lend themselves 
to a GPD interpretation, of exclusive meson electropro- 
duction). Numerous extensions of the code are possible, 
such as introducing QCD higher twists and/or higher or- 
der corrections, implementing dispersion relations, and 
so forth. 

The author has benefitted from numerous discussions 
with many colleagues from JLab and from the French 
"Nucleon GDR" group, whose list would be too long to 
enumerate but to whom the author is very thankful. Very 
special thanks are given to P. Desesquelles, M. Gargon, 
K. Griffioen, C. Mufioz Camacho and M. Vanderhaeghen. 



[1] D. Miiller, D. Robaschik, B. Geyer, F.-M. Dittes, and J. 

Horejsi, Fortschr. Phys. 42 (1994) 101. 
[2] X. Ji, Phys. Rev. Lett. 78 (1997) 610; Phys. Rev. D 55 

(1997) 7114. 
[3] A.V. Radyushkin, Phys. Lett. B 380 (1996) 417; Phys. 

Rev. D 56 (1997) 5524. 
[4] J.C. Collins, L. Frankfurt and M. Strikman, Phys. Rev. 

D 56 (1997) 2982. 
[5] K. Goeke, M. V. Polyakov and M. Vanderhaeghen, Prog. 

Part. Nucl. Phys. 47, 401 (2001). 
[6] M. Diehl, Phys. Rept. 388, 41 (2003). 
[7] A.V. Belitsky, A.V. Radyushkin, Phys. Rept. 418, 1 

(2005). 
[8] A. Airapetian et al., Phys. Rev. Lett. 87, 182001 (2001). 
[9] S. Stepanyan et al., Phys. Rev. Lett. 87, 182002 (2001). 
[10] S. Chen et al., Phys. Rev. Lett. 97, 072002 (2006). 
[11] A. Airapetian et al., arXiv:0802.2499 [hep-ex]. 
[12] A. Airapetian et al., Phys. Rev. D 75, 011103 (2007) 
[13] C. Mufioz Camacho et al., Phys. Rev. Lett. 97, 262002 

(2006). 
[14] F.-X. Girod et al., Phys. Rev. Lett. 100, 162002 (2008). 
[15] JLab approved experiment E-05-114 : "Deeply Virtual 
Compton Scattering at 6 GeV with polarized target and 
polarized beam using the CLAS detector"; spokesper- 
sons : A. Biselli, L. Elhouadrhiri, K. Joo and S. Niccolai. 
[16] JLab (conditionally) approved experiment E-08-021 ; 
"Deeply Virtual Compton Scattering at 6 GeV with 
transversely polarized target using the CLAS detector"; 
spokespersons : H. Avakian, V. Burkert, M. Guidal, R. 
Kaiser and F. Sabatie. 
[17] JLab approved experiment E-07-007 ; "Complete separa- 
tion of Deeply Virtual Photon and tt" electroproduction 
observables of unpolarized protons"; spokespersons ; C. 
Munoz Camacho, J. Roche, C. Hyde- Wright and P.-Y. 
Berlin. 
[18] JLab approved experiment E-06-003 : "Deeply Virtual 



Compton Scattering with CLAS at 6 GeV"; spokesper- 
sons : V. Burkert, L. Elhouadrhiri, M. Gargon, R. Niya- 
zov and S. Stepanyan. 

[19] JLab approved experiment E-12-06-114 : "Measurements 
of Electron- Helicity Dependent Cross Sections of Deeply 
Virtual Compton Scattering with CEBAF at 12 GeV", 
spokespersons : C. Hyde-Wright, B. Michel, C. Munoz 
Camacho and J. Roche. 

[20] JLab approved experiment E-12-06-119 : "Deeply Vir- 
tual Compton Scattering with CLAS at 11 GeV", 
spokespersons : V. Burkert, L. Elouadrhiri, M. Gargon, 
M. Holtrop, D. Ireland, K. Joo, W. Kim, F. Sabatie. 

[21] A. Belitsky, D. MuUer and A. Kirchner, Nucl. Phys. B 
629, 323 (2002) 

[22] M. Mazouz et al., Phys. Rev. Lett. 99, 242501 (2007). 

[23] M. Guidal and M. Vanderhaeghen, Phys. Rev. Lett. 90, 
012001 (2003). 

[24] A. V. Belitsky and D. Muller, Phys. Rev. Lett. 90, 
022001 (2003). 

[25] O. V. Teryaev, arXiv/0510031 hep-ph ; I. V. Anikin and 
O. V. Teryaev, Phys. Rev. D 76, 056007 (20 07). 

[26] M. Diehl and D. Y. Ivanov, larXiv:0707.03511 hep-ph. 

[27] K. Kumericki, D. Muller and K. Passek-Kumericki, Nucl. 
Phys. B 794, 244 (2008) 

[28] M. Polyakov and M. Vanderhaeghen, arXiv:0803.127T] 
[hep-ph]. 

[29] M. Vanderhaeghen, P.A.M. Guichon, M. Guidal, Phys. 
Rev. D 60, 094017 (1999). 

[30] M. Guidal, M. V. Polyakov, A. V. Radyushkin and M. 
Vanderhaeghen, Phys. Rev. D 72, 054013 (2005). 

[31] M. Guidal, M. V. Polyakov, A. V. Radyushkin and M. 
Vanderhaeghen, Phys. Rev. D 72, 054013 (2005). 

[32] F. James, MINUIT, D507, CERN (1978). 

[33] M. Polyakov and C. Weiss, Phys. Rev. D 60, 114017 
(1999). 

[34] D. S. Hwang and D. Muller, Phys. Lett. B 660 (2008) 



350. [36] F. X. Girod, PhD thesis, Universite Louis Pasteur, Stras- 

[35] N. Kivel, M. V. Polyakov and M. Vanderhaeghen, Phys. bourg, France (2006). 

Rev. D 63, 114014 (2001). 



